Decoherence of a single-ion qubit immersed in a spin-polarized atomic bath 
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We report on the immersion of a spin-qubit encoded in a single trapped ion into a spin-polarized 
neutral atom environment, which possesses both continuous (motional) and discrete (spin) degrees of 
freedom. The environment offers the possibility of a precise microscopic description, which allows us 
to understand dynamics and decoherence from first principles. We observe the spin dynamics of the 
qubit and measure the decoherence times (Ti and Tb), which are determined by the spin-exchange 
interaction as well as by an unexpectedly strong spin-nonconserving coupling mechanism. 
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A spin-1/2 system represents the most fundamental 
quantum mechanical object. Its spin-dynamics and - 
decoherence when interacting with an environment de- 
termine its potential use as a qubit and are responsible 
for a multitude of impurity effects encountered in the 
solid state. While an extensive amount of theoretical 
work on this problem exists (for reviews see [H ex- 
periments with well-controlled and adjustable environ- 
ments are scarce. The necessary sensitivity to probe 
single spins, ideally with a single-shot readout, is often 
incompatible with the ability to adjust the properties 
of the environment. Among the few examples of con- 
trolled decoherence processes are the quantum measure- 
ment process the decoherence of motional quantum 
states of trapped ions by noisy classical electric fields [U , 
and the effects of spontaneous emission 0-0| • These are 
related to spin-boson physics, where the environment is 
formed by a set of harmonic oscillator modes [lj. On 
the contrary, the decoherence of a localized spin impu- 
rity inside an environment of other spins lies at the heart 
of the so-called central-spin problem [2[. Here, the oc- 
curring decoherence mechanisms are different from spin- 
boson physics since the spectrum of the bath is discrete 
and spin-spin or spin-orbit interactions play a role. The 
central-spin model has been often applied as a simplified 
and approximate description of the interaction of semi- 
conductor quantum dots [1] and color centres in solids [§| 
with their environment. 

Here we investigate the model system of a single local- 
ized spin-1/2 coupled to an spin-polarized environment 
of tunable density. Specifically, we embed a trapped sin- 
gle Yb + ion, initially laser-cooled to Doppler tempera- 
ture, into a spin-polarized ultracold neutral bath of 87 Rb 
atoms and study the resulting decoherence of the ion's 
internal spin state. We observe an intricate decoherence 
mechanism, which is spin-nonconserving and involves the 
coupling of the orbital degrees of freedom with the spin 
degrees of freedom. We measure the longitudinal (Ti) 
and the transverse (T2) coherence times, also with respect 
to the energy separation, and identify the time scales for 
Zeeman- and hyperfine-state relaxation. 
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FIG. 1. (Color online) a Illustration of the trapped single 
ion spin coupled to a spin-polarized neutral atom cloud, b 
Level structure of the Yb + electronic ground state in a weak 
magnetic field. To implement the spin-1/2 system, we use 
either the Zeeman qubit \mj = ±1/2) in the isotope 174 Yb + 
or the magnetic field insensitive hyperfine qubit \F = 0, mF = 
0) and |1,0) in the isotope m Yb + . c The cloud of neutral 
87 Rb atoms is prepared in one of the four atomic spin states 
|F = 2,m F = 2) 0) |2 > -2)a ) H,l>o or |1,-1)„. 



The interaction between an ion and a neutral atom is, 
to leading order, via a central potential V a (r) where r 
is the internuclear separation. Asymptotically (r — > oo), 
the central potential is determined by the polarization 
interaction —C^jlr^ and is spin-independent [El. Ull. [l2| . 
For small r, collisions in the electronic singlet and triplet 
channels exhibit different potentials Vs(r) and Vr(r) giv- 
ing rise to a spin-exchange interaction, which is spin- 
conserving for the total spin of atom and ion. Langevin 
collisions in the central potential, i.e. collisions with en- 
ergies above the centrifugal barrier, occur at an energy- 
independent rate ~/l = 1 r K\jC^j fJ.n a , where fi is the re- 
duced mass and n a is the neutral atom density. In our 
experiment we have ^/i,/n a — 2.1 x 10 _15 m 3 /s and we 
use typical densities of n a = 10 18 m~ 3 . Even though the 
Langevin rate is significantly smaller than the total col- 
lision rate 7 C (see Supplementary Material), in previous 
experiments the Langevin process has dominated cold 
inelastic ion-atom collisions @, 0, EH flil - Hij ] . For spin- 
dependent processes additionally the anisotropic mag- 
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FIG. 2. (Color online) Spin-relaxation in the 174 Yb + Zeeman 
qubit. The probability pBright of the ion to occupy the bright 
| X) state after preparation in | X) (full symbols) or | \) (open 
symbols) in a bath of (a) |1, l) a or (b) |2,2) a atoms. Error 
bars denote one standard deviation uncertainty intervals re- 
sulting from approximately 3000 measurements per spin state, 
c The equilibrium spin states (corrected for detection efficien- 
cies) for all four atomic bath configurations |2, — 2) a , |1, — l)a, 
|1,1) and |2,2)„. 



netic dipole-dipole interaction and the second-order spin- 
orbit coupling can play important roles, the latter in 
particular for heavy atoms [19]. They induce coupling 
between the spin and the orbital motion and break the 
conservation of the total electronic spin. In this work we 
determine which of these mechanisms lead to spin deco- 
herence. 

First, we study spin relaxation of the Zeeman qubit of 
the isotope 174 Yb + , which represents an ideal two- level 
system (see Figure 1). We prepare the \mj = 1/2) = | f) 
or the \mj = —1/2) = | I) qubit state in an adjustable 
magnetic field inside the neutral atom cloud. During 
the interaction period, the ion undergoes binary colli- 
sions with the atoms, and we normalize interaction times 
by the Langevin time constant ti = 1/7l- We detect 
the population of the qubit state by electronic shelving 
and determining the probability of the ion being in a 
bright (p Bright) or dark (pDark) state using light scat- 
tering on the S1/2 ~ A/2 transition (see Supplemen- 
tary Material). Measurements of the T\ time are per- 
formed for varying density and neutral atom spin com- 
position. The data for the two different bath config- 
urations \F = l,rriF = l) a and |2, 2) a are shown in 
Figures 2 a and 2b, respectively. We find that the ini- 
tially polarized ion spin relaxes into a mixed steady-state 
within a few Langevin collision times. This result cannot 
be understood from spin-exchange collisions, for which, 
for example, the doubly spin-polarized combination | f) 
and 1 2, 2) a is protected from spin-changing collisions due 
to spin conservation. From our steady-state population 
data for atoms in the |2,2) a state and the measured 
T\ = (2.50 ± 0.39)£l, we determine the spin-exchange 
and spin-relaxation rates to be ji,se — (0.22 ± 0.03)/Ti 
and 7t,sfl = "fi,SR — (0.39 ± 0.02)/7\ (see Supplemen- 
tary Material). An analogous result holds for the |2, — 2) a 
state. Notably, the measured spin-relaxing collision rate 
is approximately 5 orders of magnitude higher than the 



charge-exchange collision rate 0, 13] • All measurements 
were performed at an energy splitting of the ion Zee- 
man qubit of 37.5 MHz, however the same behaviour was 
found for different Zeeman splittings in the range from 
0.39MHz to 139 MHz in a |2,2) a environment. 

The observed spin dynamics of the Yb + ion in a spin 
polarized cloud of 87 Rb atoms is inconsistent with the 
picture of dominant spin-exchange and negligible spin- 
relaxation, that has previously been reported in He + + 
Cs collisions 2l|, in neutral atom-atom collisions even 
at room temperature [22| . and in semiconductor quan- 



gases 



turn dots [23|. Moreover, investigations in both ultracold 



24] and optically pumped vapor cells [25j have re- 
vealed that even for heavy alkali-metal collisions the hi- 
erarchy -f c ~ 7se ^> 7s R is fulfilled. In such a situation, 
the dominant spin-exchange interaction leads to a steady 
state near a perfect polarization of spin and environment 
both when the environment is initially polarized but the 
spin is unpolarized and when the single spin is contin- 
uously pumped but the environment is initially unpo- 
larized. The strong spin-relaxation observed here could 
be aided by a level crossing of both the incoming sin- 
glet A 1 Y 1 + and triplet a 3 E + channels to a 3 n channel 
at short internuclear separation [26[ 27 j, which poten- 
tially provides a mechanism for spin-orbit coupling. A 
conceptually similar competition between spin-orbit in- 
duced relaxation and spin-conserving processes seems to 
play a role in the ultrafast relaxation of magnetization in 
ferromagnets [28J. 

The dependence of the steady-state spin population of 
the ion on the atomic bath configuration has been deter- 
mined for the |2, — 2) , |1, — l) tt , |1, l) a , and |2, 2) a states 
and is shown in Figure 2c. The steady-state values for the 
|1, — 1) and |1, 1) Q states can be qualitatively understood 
considering their projections into the electron singlet and 
triplet bases and taking into account that spin-exchange 
collisions leading to a transition of the atomic hyperfine 
state from the \F = l) a to the \F = 2) a manifold are 
energetically suppressed. 

Next, we study the spin relaxation of the magnetic 
field insensitive hyperfine qubit of the electronic ground 
state of 171 Yb+ (see Fi gure 1 and Supplementary Mate- 
rial). The spin detection is performed by resonant light 
scattering on the Si/ 2 \F = 1) — s- P\/2\F — 0) transition 
and measuring the probability p Bright (pDark) of the ion 
being in a bright (dark) state \F = 1) (\F = 0)). If the 
neutral atoms are prepared in the |1, — 1)„ state (see Fig- 
ure 3a), the ion relaxes towards the spin ground state, 
as expected for amplitude decoherence, and we observe 
an exponential decay with a time constant T\ of a few 
Langevin collision times. 

In contrast, if the neutral atoms are prepared in the 
|2,2) a state (Figure 3b), the steady state probability to 
measure the ion spin in the \F = 1) manifold is 0.16(1). 
We interpret this non-zero steady-state value as resulting 
from the intake of hyperfine energy from the atoms. In 
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FIG. 3. (Color online) Hyperfine spin relaxation in 171 Yb + . 
a The probability p Bright after preparation in the |1,0) (full 
symbols) or the |0, 0) (open symbols) state vs. the interaction 
time for atoms in the |1, — l) a state. Error bars denote one 
standard deviation uncertainty intervals resulting from a to- 
tal of 5000 measurements. The fit values at t — are limited 
by detection errors, b Similar data for collisions with 87 Rb 
atoms in the |2, 2} a state and a total of 19000 measurements, 
c-e Zeeman-resolved detection within the F = 1 manifold af- 
ter preparation in |1, 0) (atoms in |2, 2) a ) with 1200 measure- 
ments per Zeeman state. The measurements are performed 
by applying resonant n pulses, exchanging the population of 
the dark state |0, 0) with |1, — 1), |1,0) or |1, 1) immediately 
before the detection of the probability of the ion being in the 
dark state poark- 



order to quantitatively assess this effect, we assign the 
steady-state distribution of the ion's hyperfine spin state 

a "spin-temperature" T s by = 3exp (— ^. gT ) result- 
ing in T s = 200 mK. Here pi — Yl m P\i,m) ^ s the proba- 
bility PBright corrected for detection efficiencies (see Sup- 
plementary Material) and E HFS = hx 12.6 GHz is the 
hyperfine energy of the ion. In comparison, we estimate 
the kinetic energy Eki n from the balance of elastic and 
inelastic processes, assuming to be not dominated by mi- 
cromotion heating effects at these temperatures [29] . The 
average kinetic energy intake of the ion per Langevin 
collision due to a hyperfine flip in the neutral atom is 



SEheat = — 1J ¥ — • Here, m a is the mass of the 
atom, rrii the mass of the ion, < e < 1 describes the 
probability that an atomic hyperfine flip occurs during a 
Langevin collision, and E^ FS = h x 6.8 GHz is the in- 
ternal hyperfine energy of the atom. The average ion 
energy loss by elastic ion-atom momentum transfer is 
SE coo i = - E kin per Langevin collision. As a 



result, the steady-state average kinetic energy of the ion 
is (E km ) = e E^ FS ^±^ = ek B x 240 mK. For e near 
unity, we hence obtain \E kin ) as kgT s , which signals an 
cquidistribution of energy between kinetic and spin de- 
grees of freedom. In contrast, if the atomic bath is pre- 
pared in its hyperfine ground state \F = l) a , no energy- 
releasing hyperfine changes can occur and the temper- 
ature of ion is expected to be limited by micromotion 
heating, which is on the order of 20 mK 21, 29|. 

As an additional degree of freedom in the case of 
i7i Yb+, we study the spin transfer within the \F = 1) 
manifold of 171 Yb + when starting from |1,0). The data 
are shown in Figure 3c-e, together with the results of a 
four-level rate equation model, which involves spin trans- 
fer within the \F = 1) manifold as well as decay into the 
\F = 0) state. Most strikingly, the decay out of the ini- 
tially prepared |1, 0) state exhibits two different time con- 
stants (see Figure 3d). The fast initial decay is associated 
with the population of the |1, ±1) states, where the occu- 
pation correspondingly rises (see Figure 3c and e) . Even- 
tually, the \F = 1) states decay to their steady-state pop- 
ulations. The significant build-up of population in the 
1 1 , — 1) state and the relaxation from |1,1) — > |0, 0) are 
both forbidden by spin conservation, however, so as for 
the Zeeman qubit we observe a strong spin-nonconserving 
process. 

Quantum mechanically, an isolated spin- 1/2 can exist 
in superposition states of spin-up and spin-down that are 
not describable by classical physics. Coupling to an en- 
vironment affects the quantum mechanical superposition 
and, eventually, leads to decoherence as the information 
of the quantum correlations gets lost and the quantum 
mechanical superpositions transform into probability dis- 
tributions of occupation numbers. In the problem dis- 
cussed here, it is important to understand whether this 
decoherence results from inelastic spin-relaxation pro- 
cesses, i.e. amplitude decoherence, or whether also elastic 
(forward) scattering contributes to additional phase de- 
coherence of superposition states. In the hydrogen maser 
[30l | spin-exchange processes have been identified as a 
leading source of both line shifts and broadening. 

We characterize the influence of ion-atom collisions on 
the hyperfine clock transition in 171 Yb + by measuring 
its coherence time and line shift using the Ramsey tech- 
nique. To this end, we prepare the ion in the |0,0) 
hyperfine ground state and apply a 7r/2-pulse on the 
12.6-GHz clock transition, creating a superposition state 
(|0,0) + i|l,0))/V2. After a waiting time of 27ms we 
apply a second 7r/2-pulse, followed by a readout of the 
\F = 1) spin state. Figure 4 displays our results. In pan- 
els a-d we show sample traces of our Ramsey fringes as 
a function of the detuning of the microwave pulse from 
the clock transition with increasing density of neutral 
atoms during the interaction time. Even without neu- 
tral atoms, the visibility of the Ramsey fringes is limited 
to (55 ± 2)% owing to magnetic field fluctuations of the 
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FIG. 4. (Color online) Spin coherence of the hyperfine clock 
transition |0, 0} «-> 1 1 , 0) in m Yb + . a-d Ramsey fringes 
recorded for increasing atomic density, e Ramsey contrast 
as a function of t/tn displays an exponential decay with a 
time constant of T2 = (1.4 ± 0.2)^- f No frequency shift of 
the clock transition is observed within the experimental er- 
rors. Error bars denote one standard deviation uncertainty 
intervals. 



bias field. We find that the spin coherence decays within 
T 2 = (1.4 ± 0.2)t L (see Figure 4e). This is on the time 
scale of the population relaxation of the |1, 0) state, which 
is driven by the spin-exchange and spin-orbit interaction. 
The measurement hence identifies spin-relaxation as the 
leading mechanism of spin decoherence in our system 
with a non-detectable contribution from elastic (forward) 
scattering. Finally, we note that the frequency shift Av 
of the clock transition is below the resolution set by the 
spin relaxation rate and is Av < 4 x 10~ 11 El IFS jh for 
our densities (see Figure 4f). 

Our measurements realize a novel spin-bath with addi- 
tional continuous degrees of freedom in which we have ob- 
served an unexpectedly strong coupling between internal 
(spin) and external (orbital) degrees of freedom. Future 
experiments with a driven impurity could reveal whether 
this coupling mechanism could lead to similar quantum 
hase transitions as predicted for the spin-boson model 
1] and how the physics of the single ion in the spin- 



polarized bath relates to the physics of polarons [31H34 1 . 
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SUPPLEMENTARY MATERIAL 



Experimental setup 



We prepare up to 8 x 10 5 neutral 87 Rb atoms in the 
\F = 2,mp = 2) a hyperfine state of the electronic 
ground state at temperatures down to T w 200 nK in 
a harmonic magnetic trap of characteristic frequencies 
l] (uj x ,ujy 7 uj z ) — 2tt x (8,26, 27) Hz. When prepared 
in the \F = l,mp = — l) a state the trap frequencies 
are (lj x ,lj v ,U! z ) = 2tt x (5,10, 15) Hz. Alternatively, the 
atoms can be transferred into an optical dipole trap 
formed by two crossed laser beams at 1064 nm. Here, 
the atoms can also be trapped in the hyperfine ground 
state |1, 1) or the |2,— 2) a state. At the same location, 
we trap a single Yb + ion in a radio-frequency Paul trap 
with secular trap frequencies of uj± = 2ir x 150 kHz ra- 
dially and uJax — 2ir x 55 kHz axially @, HJ . The atomic 
density at the location of the ion in the centre of the neu- 
tral bath is determined for magnetically trapped clouds 
from atom number and temperature measured in absorp- 
tion imaging and the harmonic trap frequencies. For op- 
tically trapped clouds, the atomic density is determined 
by measuring the inelastic loss rate of the ion in its elec- 
tronically excited Z?3/2-state In both cases we esti- 
mate the uncertainty of the absolute atomic densities to 
be about 40%. In order to avoid losses of the ion due 
to inelastic collisions in the highly reactive electronically 
excited states the atoms are removed from the posi- 
tion of the ion before detection of the 174 Yb + qubit. For 
steady-state spin measurements the same atomic cloud 
has been used for up to eight consecutive measurements 
by displacing and returning the atoms by moving the po- 
sition of the optical trap beams. 



Ion-atom interaction 



The interaction between an ion and a neutral atom 
at long distances is dominated by the attractive polar- 
ization interaction potential, which is of the form of a 
central potential V(r) — — Here, C 4 = a q 2 / (Aire ) 2 
is proportional to the neutral particle polarizability ao, 
q is the charge of the ion, eo is the vacuum permittiv- 
ity, and r is the internuclear separation. The total colli- 
sion rate between atoms and ions in the central potential 
is given by 7c = n a ^(l + 7 r 2 /16)(C 4 /^) 2 / 3 (^//x) 1/6 
and includes quantum mechanical forward-scattering |5| . 
Here, E is the collision energy. In our experiment we have 
lL /n a = 2.1 x 10~ 15 m 3 /s and j c /n a = 2.5 x 10~ 14 m 3 /s 
for a collision energy of E — ks x 100 mK. 
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Ion qubit preparation and detection 

Ultraviolet light at 370 nm, resonant with S1/2 P\ii 
transition, is used for laser cooling and spin state prepa- 
ration, while laser light at 935 nm clears out the 2 D 3 / 2 
state. State preparation in either Zeeman state of the 
S1/2 ground state of the 174 Yb + isotope is achieved by 
optical pumping with er + — or (7~— polarized light on the 
<Si/2 — ^ P1/2 transition. State detection is performed 
by frequency-selective electron shelving in the F7 / 2 state 
via the -D5/2 state using a laser at 411 nm and subse- 
quent resonant light scattering on the Si /2 — > Pi /2 tran- 
sition Imperfections in the detection scheme limit 
the direct correspondence of the detection of a bright 
or dark quantum state, as determined from the statis- 
tics of photon events, and the occupation of a specific 
upper or lower spin quantum state. The probability 
of spin occupation p^ is determined by p^ = (poark ~ 
VDark,t)/(ilDark,i - VDark,f) from the detected distribu- 
tion of dark states pDark and the independently deter- 
mined detection errors (1 — t]Dark,i) and r]Dark,t- F° r 
measurements with 87 Rb \F — 2) a , the detection errors 
are (1 - r) Darkti ) = 0.19 ± 0.01 and r) Dar k,t = 0.03 ± 0.01, 
and for measurements involving the |_F = l) a manifold we 
find (l-7]Dark,i.) = 0.10±0.01 and r, Dar k,t = 0.00 + 0.01. 
The larger detection errors of the Zeeman qubit when the 
atoms are in the \F — 2) a state result from the increased 
ion kinetic energy due to release of atomic hyperfine en- 
ergy- 

For the isotope 171 Yb + , state preparation in the |0, 0) 
state is achieved by a short (250 ps) laser pulse reso- 
nant with the Si/2\F = 1) — » P1/2I-F = 1) transition 
Microwave radiation resonant with the 12.6-GHz 
hyperfine transitions is used for coherent qubit manip- 
ulation on the clock transition, for Rabi 7r-flips on the 
|0,0) -» 1 1 , — 1 ) and |0,0) -> |1,1) transitions as well 
as for repumping during Doppler cooling. The mag- 
netic bias field strength depends on the atom trap con- 
figuration and ranges from 5.64 G in the optical trap to 
7.56 G in the magnetic trap of the 87 Rb|l,— l) a atoms. 
The spin detection probes fluorescence of the ion on the 
S1/2 \F = 1) — > P1/2 \F = 0) transition and we use a 



standard statistical analysis of photon detection events 
for state discrimination |7J. The fluorescence detection 
can distinguish the states |0, 0) and \F — 1) and the 
detection errors are (1 — rjD a rk,i) = 0-02 ± 0.01 for the 
hyperfine ground state \F = 0) and r]Dark,t = 0.07±0.03 
for the \F = 1) hyperfine states. 

Determination of spin-exchange and spin-relaxation 

rates 

The spin relaxation behaviour is well described by 
solutions of the general two-level rate equation model 
p t (t) = p tfi e- t / Tl + pt,«>(l ~ e~ t/Tl ) with the lon- 
gitudinal coherence time T\ = 1/(7-1- + 74-) an d the 
equilibrium spin state Pf j0 o = 7t/(7t + 74. )■ The de- 
cay rates for the spin states 7^ = 7^,sb + 7t,Sfl and 
74 = 74, SB + li,SR are decomposed into a spin-exchange 
part depending on the atomic bath spin and a spin- 
relaxation part including spin-orbit coupling. For col- 
lision with the maximally polarized bath state |2,2) Q at 
energies much larger than the Zeeman splittings, the re- 
lations 7f.s£ = and 7^s_r = J^,sr hold true and re- 
sult in j itSE = (2p TiCO - 1)/Ti ='(0.22 ± 0.03)/^ and 
7t,sfl = H,sr = (1 -Pt.ooVT 1 ! = (0.39 ± 0.02)/T x . 
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Ion isotope 


Atom state 


Ti/tz, 


Pt.oo 




Tilth 


171 Yb+ 




1.73 ±0.17 




0.000 + 0.005 




171 Yb+ 


|2,2)„ 


3.39 ±0.16 




0.163 ±0.013 


1.4 ±0.2 


174 Yb+ 


|2,2) a 


2.50 ±0.39 


0.609 ±0.015 






174 Yb+ 


|2,-2)„ 




0.423 ± 0.026 






174 Yb+ 




1.60 ±0.24 


0.563 ±0.017 






174 Yb+ 


|l,-l)a 




0.457 ±0.021 







TABLE I. Summary of the decoherence times (Ti and T2) and steady-state spin distributions of all spin-bath combinations. In 
addition to the statistical errors quoted, all times have systematic uncertainties of 40% due to the determination of the absolute 
atomic densities at the location of the ion. All values are compensated for detection efficiencies. 



